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ABSTRACT 

We present new analytic calculations of the coupling between ultraviolet resonance photons and the pop¬ 
ulation of the hyperfine states in the light elements (H, D, 3 He + ) which include several previously neglected 
physical processes. Among these are the backreaction of resonant scattering on the pumping radiation, the 
scattering of Ly/3 photons and the effect of local departure from pure Hubble flow. The application of the new 
treatment to the redshifted hydrogen 21 and deuterium 92 cm lines from the high-redshift universe results in 
an amplitude correction of up to an order of magnitude. We further show that the standard assumption that 
ultraviolet pumping drives the spin temperature towards the kinetic temperature does not hold for deuterium, 
whose spin temperature is generally negative. 

Subject headings: cosmic microwave background - cosmology: theory - diffuse radiation - intergalactic 
medium - radio lines: general 


1. INTRODUCTION 

The hyperfine transitions in hydrogen, deuterium and 
helium-3 provide an opportunity to study the universe during 
its “dark ages” and the epoch of reionization that followed. 
From the radio continuum produced by redshifted lines of 
these elements, which can be observed in absorption or in 
emission against the cosmic microwave background (CMB), 
information about the thermal evolution, density perturbation 
spectrum and ionization history of the universe can be ob¬ 
tained. However, the correct interpretation of the measured 
signal requires knowledge of the processes which decouple 
the spin temperature of the atoms, 7' s , from the CMB tem¬ 
perature, 7cmb- Of these, one of the most important during 
the epoch of reionization (6 < z < 20) is ultraviolet pump¬ 
ing. Resonant photons, such as Lya, absorbed by the atom 
in the ground state may, after the decay of the excited state, 
leave the electron in a different hyperfine level. This process 
couples the spin temperature to the color temperature of the 
resonant photons, T a (Field 19581) 


j _ Tcmb +y a T a +y c Tk 
1 +y a +yc 


( 1 ) 


where 7 k is the kinetic temperature, y c is a constant propor¬ 
tional to the collisional excitation rate and y a is given by 


ya = 


P\qT* 
A i 0 T a 


( 2 ) 


where I, = Invhyp/k, Uhy P is the frequency of the hyperfine tran¬ 
sition, A io is the spontaneous emission coefficient of the hy¬ 
perfine transition and P;o is the radiative deexcitation rate for 
the upper hyperfine level due to Lya photons. 

As we show in this paper, the correct application of the 
above expressions requires us to take into account several 
physical processes which have previously been neglected. 
The first of these is the effect of the hyperfine splitting on the 
photon color temperature and intensity. In their pioneering pa¬ 
pers, Wouthuysen (1952) and Field (1958, 1959) argued that 
in the limit of large optical depth, T a must approach the ki¬ 
netic temperature of the atoms, to which the photons are cou¬ 
pled by recoil. As a result, since the Lya resonance scattering 


depth for hydrogen atoms in the neutral intergalactic medium 
(IGM) is extremely large, it has been customary ever since 
then to replace T a with 7 k in eqs. 1 and 2 above. That argu¬ 
ment, however, did not take into account the effect of the hy¬ 
perfine splitting on the change of photon energy during scat¬ 
tering. Just as the spin temperature of the atoms is affected by 
the spectrum of the UV photons in the resonance line, so the 
UV spectrum around the resonance is in turn affected by the 
spin temperature. This backreaction shifts T r , towards 7 cmb, 
thus reducing the coupling between the kinetic and spin tem¬ 
peratures. In addition, we consider the effect of scatterings on 
the radiation intensity near the resonance. Chen & Miralda- 
Escude (2004) showed that, in a homogeneously expanding 
cosmological medium, the intensity of the continuum radia¬ 
tion drops around the resonance, thus reducing the value of 
y a . Here we derive an analytical formula for the amplitude 
of this drop. We find significant discrepancies between our 
analytical calculations and some of their numerical results. 
Furthermore, we extend their analysis to consider departures 
from pure Hubble flow, associated with cosmic structure for¬ 
mation, such as a contracting medium which is appropriate 
for overdense regions. Finally, we also consider scattering by 
resonant photons other than Lya. Lya photons have the high¬ 
est scattering cross-section and, unlike higher resonances like 
Ly j3, they are not destroyed by cascade after just a few scat¬ 
terings. For these reasons, the higher resonances are generally 
assumed to be unimportant. While correct for hydrogen, this 
assumption proves to be a very gross error in the case of deu¬ 
terium, for which the role of Ly/3 photons is crucial. 

The corrections derived here to the hydrogen spin tem¬ 
perature which results from Lya pumping in the early 
universe can have a significant effect on predictions of the 
21-cm brightness temperature and its fluctuations. As we 
shall show their importance depends on the intensity of 
the UV background and the temperature of the IGM. An 
approximation which is often made in recent theorecti- 
cal predictions of the 21-cm background from the epoch 
reionization involves the assumption that the UV pumping 
background is high enough to make T s equal 7)., while 
some other process (e.g. X-ray background) heats the IGM 
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21-cm brightness temperature becomes independent of the 
actual value of T s , so the calculations are greatly simplified. 
However, the period either before the UV background reaches 
such high intensity or the neutral IGM is heated is not well 
represented by this limiting approximation. This includes the 
epoch when the 21-cm signal is in absorption and potentially 
at its strongest. That is the case in which the results presented 
here will be important. 

While this paper was in preparation, a preprint of Hirata 
(2006) was posted in which some of the above effects were 
calculated numerically. We find an excellent agreement be¬ 
tween our analytic solutions and his numerical results. 

In §2, we estimate the effect of the hyperfine transition on 
the color temperature of Lya photons. In §3, we calculate the 
change of intensity of Lyo photons due to their scattering by 
atoms. In §4, we estimate the effect of these processes on the 
radio emission signals from D and 3 He + . In §5, we discuss the 
implications of our results for the interpretation of hydrogen 
and deuterium radio signals from the epoch of reionization 
and the end of the cosmic dark ages. 


2. THE COUPLING BETWEEN THE UV AND RADIO SPECTRA 


Consider a UV resonance transition which is followed by 
a single-photon decay to either hyperfine level of the ground 
state. Upon such scattering by an atom, the energy of the in¬ 
coming UV photon is changed due either to atomic recoil or 
to the hyperfine transition. In most cases, both the kinetic and 
the spin temperature of the atoms is below the energy of the 
photons ( hv/k ~ 10 5 K), so on average the photons will lose 
energy in scatterings. The actual energy change, however, de¬ 
pends on the exact frequency of the photon. Assuming that the 
energy distribution of atoms is Maxwellian, we found that, in 
each scattering, the average energy loss of the photons is (see 
Appendix for details) 


A Ei = 


(Ivy ) 2 


me* 


, kT k a'(v)\ 
li u(v) ) 


due to the atomic recoil and 


A E,= 


b(hv kw ) 2 

2kE 


1- 


kT s a'(v) 
h cr(i/) 


(3) 

(4) 


due to the change of the hyperfine state, where <r(v) is the 
scattering cross-section, a'(v) = da/dv, is the frequency 
of the hyperfine transition, b is the probability that a single 
resonance scattering results in a spin-flip and m is the atomic 
mass. Similarly, the total energy loss in the two processes 
may be written as 

AE = AE l+ AE 2 = (5) 

where A = (h 2 v 2 /mc 2 ) + (bli 2 i'l yp /2kT s ) and B = 
(T k h 2 u 2 /mc 2 ) + (bh 2 vl w /2k). 

Consider the effect on the incoming radiation spectrum of 
repeated scatterings like those described above. This is rele¬ 
vant whenever the incoming photon is scattered without the 
possibility of photon destruction by cascade. For the radia¬ 
tion spectrum, increasing A and B by the same factor produces 
the same effect as an increase in the optical depth. Thus the 


color temperature of the resonant photons, which approaches 
a constant in the limit of high optical depth, can depend only 
on the ratio of A and B. Considering either the recoil or the 
hyperfine transition alone would give us either B/A = T k or 
B/A = T S , respectively, so obviously B/A = T c . For hydrogen 
Lya, b = 2/9, and the color temperature is given by 


(i+2.5r k )7; 

1+2.57; 


(6) 


where all temperatures are in units of degrees Kelvin. Pre¬ 
dictably the value of T a is between 7) and 7cmb- Solving 
equations 0 and © then gives the correct spin temperature 


T s = 


7cMB + [>’a,e// + >'c]7k 


1 + y (»,«//+yc 




P10T * 


MoT k 


1+0.47T 


(7) 

( 8 ) 


The the result in eq. |7]is identical to the commonly used ap¬ 
proximation for 7; described in §1, except that the value of 
the coupling constant y a is reduced by a factor (1 +0.47/7 1 ). 
In the next section we will show that the radiative deexcitation 
rate, / J m, is also affected by the changing spectrum. 


3. THE INTENSITY OF LYa PHOTONS 
3.1. Expanding medium 

The Lya pumping photons can be divided into two groups 
according to their origin. The first group (“continuum pho¬ 
tons”) consists of continuum photons originally emitted in the 
frequency range between the Lya and Ly/3 transitions, which, 
due to the cosmological redshift, are eventually shifted into 
the Lya resonance. The second group (“injected photons”) 
owes its origin to photons originally emitted between Ly 7 and 
Ly-limit frequencies. When these photons redshift into one 
of the Lyman series resonances, they get absorbed by atoms, 
which are thereby excited to n > 3 state. Since the subsequent 
cascade mostly likely goes directly to the ground state, most 
photons are reemitted close to their original frequency (with a 
small change due to a Doppler shift) and are soon reabsorbed. 
However, after just a few scatterings most of the original pho¬ 
tons are destroyed by splitting into several photons. Roughly 
two thirds of all cascades from n > 3 states and all off the cas¬ 
cades from n = 3 state, go through the intermediate 2s state, 
which decays by two-photon emission. The rest pass through 
the 2 p state, whose decay produces a Lya photon. Unlike the 
“continuum photons”, which reach the Lya resonance from 
the blue wing, the“injected photons” are produced directly in 
its center. 

When the Lya optical depth is high, the spectrum of the 
photons near the Lya resonance can be found using the 
Fokker-Planck approximation. In the comoving frame of the 
gas, the intensity J( v) varies with frequency in the neighbor¬ 
hood of the resonance at v a according to the following (see 
Appendix B for derivation) 

4>(x)J'(x) + 2[r/<p(x)+'y]J(x) = 27 / 0(1 —k a Q(x)), (9) 

where x = (v/v a — \)/(2kT k /mc 1 ) 1 l 2 is the frequency distance 
from the line center in dimensionless units, cj>(x) is the nor¬ 
malized absorption cross-section, Jo is the intensity on the far 
red side of the resonance (i.e., where it is unaffected by scat¬ 
tering) of either the continuum or injected photons, and k a 
equals zero for continuum or unity for injected photons. 

If hyperfine transitions are neglected, the recoil parameter, 
p , equals hv a /(2kT k mc 2 ) 1 ! 2 iBasko 19811) and 7 -1 = tgp is the 
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Gunn-Peterson optical depth to Lyo resonance scattering. In¬ 
cluding the effect of the hyperfine transition gives a correction 
of order (1 + 0.4/ ) to the values of 7 and 1 / (see Appendix 

B). Generally, however, this correction results in just a few per 
cent change in the radiation intensity (see Fig. |TJ. Neglecting 
terms of order 7 , the solution of equation © can be written 
as 1 

J(x ) =./( 0)e~~^ 2r,x , for k Q = 1 and x > 0, 


■ Joe 3 ° 2r,x / 27T7A l e 3 » +2r,z z 2 dz, otherwise, (10) 


where a = A2i(2kT^/mc 2 )~ 1 ^ 2 /4nv a and A 21 is the Einstein 
spontaneous emission coefficient of the Lyo transition. When 
the optical depth is low or the temperature is high, so that ther¬ 
mal width of the line becomes similar to the width of absorp¬ 
tion feature, the approximations we used in deriving equations 
®-G3 break down and higher order terms become more im¬ 
portant (iRvbicki 20061 iFurlanetto & Pritchard 2005) . How¬ 
ever, under such conditions the effect of Lyo- scattering on the 
radiation intensity is in any case negligible and the equation 
m still yields the correct answer. 

The kinetic temperature of the gas, '4, is coupled to 7 cmb 
by inverse Compton scatterings at redshifts z > 150. There¬ 
after, until the IGM is reheated during the epoch of reion¬ 
ization, the gas temperature falls adiabatic ally, according to 
7 k ss 0.02(1 + z) 2 . Prior to reionization, the optical depth of 
the neutral hydrogen evolves according to 


t gp = 7-10 5 


( ^bh 100 

^ 0.03 



l+z \ 
10 ) 


1/2 


— 

H 


-1 




( 11 ) 


where H' and 5 are the local values of the Hubble constant (to 
account for departure from unperturbed Hubble expansion) 
and the overdensity, respectively. Figure Q shows an exam¬ 
ple of the spectrum profile, J(x), for conditions corresponding 
to unheated and unperturbed neutral hydrogen at z ~ 12. Pre¬ 
dictably, near the center of the line, where frequent scatterings 
with atoms make photons lose energy faster, J(x) shows a dip, 
whose amplitude is the same for the continuum and injected 
photons 


Jo y/3 

where £ = {\6rf‘a/9it')) l l 2 (for 


■ +1 Fi ( 1 ; 


1 2 c 2 


3 , 3 , 4 ,, ( 12 ) 

hydrogen ( = 
4.6 ■ HrV l/ 2 T k -'), iF 2 is a hypergeometric function 
and J 1 and J_ are the Bessel functions of the first kind. 


3 3 

Conveniently J(0)/Jo has the same value for continuum and 
injected photons. 

Typically, (j is well below unity, which allows us to approx¬ 
imate the eq. II 2 lhv 


j_m = e - 1 .69e'\ (13) 

7o 

Since most of the scatterings happen around x = 0, the overall 

scattering rate and y a fall by the same factor, e ~'- 69 C 2/3 . Hi- 
rata (2006) found this factor numerically, separately using a 
Fokker-Planck approximation and Monte Carlo simulations. 
Figure 0 shows excellent agreement between our equations 
ll2land[T31 and Hirata’s numerical fit. 


1 It has been recently pointed to us that a similar solution for an expanding 
medium was found by Grachev(1989). 



Fig. 1.— The spectrum profile around hydrogen I.yo resonance in an 
expanding region with tgp = 10 6 and 7 - 3 K and in a contracting region 
with tgp = 3 ■ 10 7 and 7 = 10 K. The solid and dashed curves are for T s = 7f. 
for continuum and injected photons respectively. The dotted and dashed- 
dotted lines are for T s = 7j y!B = 35 K, for continuum and injected photons 
respectively. 



Fig. 2.— Hydrogen optical depth tgp for Lyo scattering at z = 12 vs 
overdensity, S = (p/ < p >)— 1. 


3.2. Contracting medium 

Equation 0 may after a slight change be used for a con¬ 
tracting, rather than expanding, medium, according to 

(f>(x)J'(x) + 2[r](f>(x) — 'y]J(x)= 2'yJo(k a O(—x) - 1). (14) 


In a contracting gas, the energy loss from scatterings slows 
down the blueshifting of the photons. Therefore instead of 
falling J(x) actually rises above Jo at the center, according to 



+ 1 Cl l; 


1 2 _ C 2 

3 ’ 3’ ~4 


(15) 


where Bi is the Airy function. Once again, this result applies 
to both continuum and injected photons. 


4. THE RARE ISOTOPES 


















4 



Fig. 3.— /(0)/7o at z = 12 for <5 = 0. The solid, dashed and dashed-dotted 
lines are obtained using respectively, eqs. I 12 l and f 01 and the numerical fit of 
Hirata (2005) for T s = Tqmb (upper lines) or T s = 7 \ (lower lines, where solid 
and dashed-dotted lines in this case are indistinguishable). 



Fig. 4.—The spectrum profile around deuterium Ly (3 resonance for 7 = 0.1 
(solid line) and 0.2 (dashed line). 


4.1. Deuterium 

When absorption of a Ly/3 photon excites a hydrogenic 
atom to the 3p state, the decay proceeds either directly to the 
Is state, in which case another Ly/3 photon is emitted, or via 
2s state, in which case the original Ly/3 photon is split into 
three photons. At redshifts z < 20, the optical depth for deu¬ 
terium Ly/3 photons is below 10. Since only 1 in ~ 8 of the ab¬ 
sorbed photons is destroyed by three-photon reemission, the 
ratio of Ly/3 to Lya scatterings does not drop as drastically 
for deuterium as for hydrogen, whose Ly/3 optical depth is 
much higher. Still, this ratio is rather low (~ 0.1 if the UV 
continuum spectrum is nearly flat between Lya and Ly 7 ), so 
Ly/3 photons would not have been important for deuterium ei¬ 
ther, if their color temperature were determined by the same 
process as for Lya . Unlike Lya photons, however, whose 
energy is changed only slightly by the scatterings and whose 
spectrum near the resonance is relatively flat, the spectrum 
around the Ly/3 resonance is determined predominantly by 
photon destruction. As photons are redshifted across the reso¬ 
nance, more and more of them are destroyed, thus making the 
red wing of the resonance much weaker than the blue. 

Since the Ly/3 optical depth of deuterium is not very high, 
the Fokker-Planck approximation can no longer be used, and 
the integral equation must be solved instead. For an isotrop¬ 
ically expanding medium, the transfer equation for radiation 
in a steady state takes the form 

0 (x)J(x) — ^J'(x) = J R(x,x )J(x )dx , (16) 

where R( x.x') is the redistribution function 
feybicki & Dell’Antonio 19941 . For the photons which 
survive the scattering (i.e., for those which are not destroyed 
by cascade) R(x,x') = Erfc(Max[|x , [, |x|])/2 iUnno 1952i) . 
Taking into account the fraction of the photons that 
are destroyed, f\, gives the full redistribution function 
R(x,x') = (1 -/d)Erfc(Max[|x , |, |x|])/2. For moderate op¬ 
tical depth the normalized absorption cross-section can be 
approximated as <f>(x) = e~ x I \pn. Using this we obtain 

e —=J{x)- 0 j'(x) = 

V't 


1 

2 


(1 - /d)Erfc(Max[|x'|, \x\])J(x')dx'. 


(17) 


Solving equation numerically (see fig. 0 , we found 
that, around the Ly/3 resonance, the color temperature is Tg ~ 

(-0.2T k 1/2 ) K. 

If the UV radiation field is sufficiently strong that we can ig¬ 
nore both the CMB photons and collisional coupling, the spin 
temperature is determined by the relative numbers of scat¬ 
terings of Lya and Ly/3 photons with deuterium atoms (the 
combined effect of Ly 7 and higher resonances is generally 
a few times smaller than that of Ly/3). Taking into account 
the higher scattering cross-section of Lya photons and the de¬ 
crease of Ly/3 intensity at the resonance, we obtain 


T s ~-IK T^ /2 n a g, (18) 

where n a g is the ratio of Lya to Ly/3 photons. For hot radi¬ 
ation sources, which produce most of the UV photons, n a g 
is close to the ratio of the bandwidths producing the photons 
(v Ly g - i^Lya)/(^Ly-y ~ ^Ly/ 3 ) ~ 3. Thus deuterium can appear in 
emission with respect to the CMB radiation even if its kinetic 
temperature is below the CMB temperature. 


4.2. 3 He 4 

The case of 3 He 4 emission is a relatively simple one. Since 
the corrections to radiation intensity and color temperature 
that were estimated in §2 and §3 are important only at low 
kinetic temperatures, we can neglect them for ionized he¬ 
lium which becomes abundant only at T > 10 3 K. Further¬ 
more, unlike deuterium, 3 He + resonances are to the red of 
4 He + . Thus practically all resonant photons except Lya are 
prevented from reaching 3 He 4 , being instead absorbed and de¬ 
stroyed by 4 He + ions first. 


5. DISCUSSION 

We have shown that the accurate modeling of the photon 
spectrum around the hydrogen Lya resonance results in a sig¬ 
nificant correction to the radiative pumping efficiency. For 
the cosmological parameters given by the latest WMAP re¬ 
sults, flboho = 0.03 and O,„o = 0.25, dSnergel et al. 20061) and 
unperturbed Hubble expansion, T s is correctly given by eq. |7] 












5 



Fig. 5.— Hydrogen spin temperature at z— 12 for different radiation inten¬ 
sities, y a o7k = 1» 10, 10 2 , 10 3 K (lower curves correspond to higher values 
of y a o). The dashed and the dashed-dotted lines are calculated, respectively, 
with and without the correction to the coupling constant. The solid lines are 
Tcmb an d ?k. 



Fig. 6 .— The correction factor to the gas absorption signal for differ¬ 
ent radiation intensities y a o^k = 1 (solid), 10 (dashed), 10 2 (dotted) and 10 3 
(dashed-dotted line). 


if y a ,eff is given by 


ya, e ff _ „-0.37(l+zW z 2t 2/3 

- — e K 



(19) 


where y a q is the previously used coupling constant, given by 


167r 2 r^e 2 /i 2 7o 

21Ai 0 T k m e c 


( 20 ) 


where / 12 = 0.416 is the oscillator strength of the Lyo transi¬ 
tion and Jq is the intensity at Lyo: resonance, when the back- 
reaction caused Lyo scattering by is neglected . 

Prior to reionization the mean differential brightness tem¬ 
perature of the hydrogen 21 cm signal is given by 


T s ~ 7c mb \ 

“iH x 

n b hm \ ( On V 1/2 f 1 + z\ 1/2 

0.03 ) \0.25 ) V 10 y 


ST b = 0.03 K 


( 21 ) 


From the above equation it follows that the strongest signal 
is obtained when the spin temperature is well below 7c MB , 
in which case T b oc T~ l . When radiation sources begin to 
form, 7 S indeed starts to drop as it decouples from 7cmb 
and approaches 7k instead. However, if X-ray photons are 
present, with an intensity which increases simultaneously 
with those in the UV, these would at the same time heat the 
gas, so that, eventually, 7 k and T s become greater than 7c MB . 
Chen & Miralda-Escude (2004) estimated that in this case the 
strongest signal f<57b ~ -20 mK) arises when the gas is still in 
absorption and T s is not yet fully decoupled from 7 cmb, i.e., 
when the signal is still very sensitive to the changes in the 
coupling strength, y a . 

A further complication arises from the small-scale density 
fluctuations. Even though the resolution of planned future ex¬ 
periments is currently limited to objects of at least ~ 1 Mpc 
size, which are still mostly in the linear regime at the rele¬ 
vant epoch, the spectrum around the resonance is determined 
by conditions on much smaller (< 1 kpc) scale, where non¬ 
linear departures from Hubble flow are important. To illus¬ 
trate the signal dependence on the small scale-structure we 


have plotted the spin temperature of the adiabatically evolv¬ 
ing hydrogen gas at z = 12 at various overdensities (Figs. |2j 
|6j(. In general, there is no one-to-one relationship between 
the density and the expansion rate, so to make this figure we 
assumed that all over/underdensities evolve as spherical “top- 
hat” density perturbations (which explains the break at 6 = 4.6 
where the gas switches from expansion to contraction). This 
yields a unique dependence of tgp on local overdensity for 
a given redshift (Fig. Q The results are shown in Figure 0 
for different values of y a fiTk (be. different intensities), which 
can be related to the ratio of photons in the UV background 
in the frequency range between Lya and Ly-limit per atom, 
N a , according to y a fiTk ~ 0.2ALU + z) 3 . The spin temper¬ 
ature changes significantly when the corrections to the cou¬ 
pling constant (eq. na are taken into account, except when 
the radiation intensity is very high. Figure [6] shows how the 
local contributions to the overall 21 cm signal change after in¬ 
cluding this correction. The net effect of departures from uni¬ 
form Hubble expansion on the mean 21 cm signal depends on 
the weighted contributions from gas at different overdensities, 
and, since this dependence is non-linear, the contributions of 
over- and under-dense regions do not cancel each other out. 

For hyperfine transition of deuterium, we found that a dras¬ 
tic difference is made by inclusion of Ly/J photons in the anal¬ 
ysis, which practically do not affect hydrogen. Thus the ratio 
between the deuterium and hydrogen radio signals is not con¬ 
stant, but depends sensitively on the UV spectrum. Unfor¬ 
tunately, this eliminates the possibility of the precision radio 
measurementsof D/H ratio suggested by Sigurdson & Furlan- 
etto (ISi gurdson & Furlanetto 20(H) . except for the period be¬ 
fore the first UV sources turn on (z. > 40). 
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APPENDIX 

ENERGY EXCHANGE BETWEEN PHOTONS AND ATOMS 

We consider the scattering of the photon of frequency v from the atom approaching with velocity v. Let us first ignore the 
hyperfine splitting. After the scattering the photon energy changes by 


v his 

hAv = hv( -y)(l -cosai), 

c me 1 


(Al) 


where a\ is the angle between direction of the photon before and after the scattering. In the rest frame of the atom, the photon 
frequency before scattering is i/(l -vcosct 2 /c), where ao is the angle between the direction of the photon and of the atom, so the 
average loss in energy averaged over all directions and velocities is 


A£i = hv 


f / tv 

/ (f>(v( 1 — vcos a .2 /c))(-^)( 1 - cos a i )P(a i, ao , v)da\da. 2 dv, 

J c me 


(A2) 


where (j)(v) is the normalized scatteri ng c ross-section and P(a \. ao , v) is the distribution function. For isotropic radiation field 
and Maxwellian velocities integrating|A2]gives 

A (hv) 2 kT k 6'(v) 

AEi = ^4(l--^^), (A3) 

me- n <p(v) 

where we have used <fi(v(\ — v/c)) « (f>(v) + <j)'(v)vv/c. 

Now lets consider the effect of the hyperfine splitting. Let b be the probability that the electron excited from lower/upper 
hyperfine level decays after photon reemission to a different hyperfine level. Further let a\(v) and <j u (v) be the photon absorption 
cross-section for the atoms at respectively low and upper hyperfine state. Then the average energy loss of a photon in a single 
scattering caused by change of electrons hyperfine state is 


- <r u 

A Ei = bhvh yp ( : —) 


(A4) 


<j\e h ^yp/ kT ‘ + cr u 

Since the frequency of the hyperfine transition, v kyp , is very small, we can approximate cr u (v) = a\(u + ^hyp) ~ cri(i/)+ v kyv (j[(v), 
<j[{v)/<j\{v) « (j>'(v)/(f>(v) and e Aiy i>yp / kT ‘ « 1 +/ir'h yp /£7’ s . Thus lA4l becomes 


A E 2 - 


b(hu h yp )- 
2kP 


h ) ' 


(A5) 


RESONANCE PROFILE 

When resonant scattering and Hubble expansion are the only important processes, the spectrum evolves according to 


dJ(v) dJ(v) 


Hv+nc / [a(iy— Av)J{y — Av) — a(v)J(v)]P(/S.i')d(Av) 1 


(Bl) 


dt dv 

where P( Av) is the probability of photons frequenc y to change by Av in a single scattering and n is the scattering atoms density. 
Neglecting terms of order (Azz) 3 we can rewrite ea. IB li as 


d.Hv) dJ(v) d 2 {J(v)(j(v)) (Av ) 2 d(J(v)(j(v)) 

-Hu+nc[ -—^-< —-— >---< Av >], 


(B2) 


dt dv **" '“ L dv 2 2 dv 

where < Av > and < (Av) 2 > are the expectation values for Av and (Av) 2 . Based on Appendix A we can show that close to 
the resonance 


< Av > = 


< (Av) 2 > - 


hv 2 


{ kT k <t)'(v)\ | bhv 2 yp 

h (f>(v) J 2kT s 


] kT s c/>'(v)\ 

h (j)(v) J ’ 


2 v 2 kT k 


+ ^%p + 


hv 2 


(B3) 

(B4) 
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The third term on the right side of the eq. IB4l can be neglected, since for 7 k > 1 K, it is at least ~ 1 () ' times smaller than the first 
term. Since the time photons spend close to the resonance is much shorter than a Hubble time, near the resonance the spectrum 
evolution can be well approximated by a steady state, i.e., dJ(v)/dt = 0. Thus we can rewrite the ea. lB2l as 


d_ 

dv 


J{v)Hv + 




, bh % P \ 

me 2 2kT s J 


<j{v)J(v) + 




a(v) 


dJ(v) 

dv 


= 0 , 


(B5) 


Integrating the above equation and expressing it in dimensionless units gives 

<t)(x)J'(x) + 2[rj(/)(x) + j]J(x) = A , (B 6 ) 

wherex= (u/u a -l)/(2kT^/mc 2 ) 1 / 2 , r] = (l+w/T s )(l+w/T^)~ l hu a /(2kT] li mc 2 ) 1 ^ 2 , w = bvl yp mc 2 /2v 2 k, 7= rQp(l +w/T s )~ i and 
tgp is the Gunn-Peterson optical depth to Lyo resonance scattering. When photons are far away from the center of the resonance 
(i.e., |x| 3 > 1) their intensity is unaffected by scatterings. Thus, since J'(x) and <p(x) go to zero in the limit |x| 3> 1, A /27 must be 
the UV intensity far away from the center of the resonance. 

As ex pected, if the probability of the spin-flip, b, is artificially set to be zero (i.e., the effect of hyperfine transition is neglected), 
eq. IB 6 l becomes identical to the one derived by Grachev (1989), based upon Chugai(1987), and later rediscovered by Chen & 
Miralda-Escude (2004), based upon Rybicki & Dell’Antonio (1994). 






